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Abstract 


In these last few decades, control charts have received a growing interest because of the 
important role they play by improving the quality of the products and services in industrial 
and non-industrial environments. Most of the existing control charts are based on the 
assumption of certainty and accuracy. However, in real-life applications, such as weather 
forecasting and stock prices, operators are not always certain about the accuracy of an 
observed data. To efficiently monitor such processes, this paper proposes a new cumulative 
sum (CUSUM) X chart under the assumption of uncertainty using the neutrosophic statistic 
(NS). The performance of the new chart is investigated in terms of the neutrosophic run 
length properties using the Monte Carlo simulations approach. The efficiency of the pro- 
posed neutrosophic CUSUM (NCUSUM) X chart is also compared to the one of the classical 
CUSUM X chart. It is observed that the NCUSUM X chart has very interesting properties 
compared to the classical CUSUM X chart. The application and implementation of the NCU- 
SUM X chart are provided using simulated, petroleum and meteorological data. 


Introduction 


Control charts play an essential role in monitoring processes in the production and 
manufacturing sectors. Shewhart-type control charts are the most popular charts because of its 
simplicity and attractive sensitivity towards large shifts in the process parameters. These charts 
are widely applied in many fields and operations (such as manufacturing processes, health care 
monitoring, stock exchange, credit card, financial fraud detection, weather updates, and inter- 
net traffic flow, etc.). However, Shewhart-type charts are inefficient in detecting small to mod- 
erate process shifts. Due to this drawback, researchers have developed new control charts 
based on classical and adaptive sampling designs to detect small-to-moderate process shifts as 
quickly as possible, including large shifts. Control charts are classified into memory-type (e.g. 
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cumulative sum (CUSUM), exponentially weighted moving average (ЕММА), etc.) and mem- 
oryless (e.g. Shewhart-type) control charts; see for instance [1]. The latter uses only the latest 
information of the process to compute the charting statistics, which makes it less sensitive to 
detect small-to-moderate shifts, but enhances its ability to detect large shifts in the process. 
The first (i.e. memory-type control chart) uses both current and past information to compute 
the charting statistics, which makes it more sensitive in detecting small-to-moderate shifts and 
less sensitive in detecting large shifts in the process. The CUSUM and ЕММА charts are the 
most popular memory-type control charts. To enhance the aforementioned charts, various 
advanced techniques are used including the addition of run rules, the use of variable sampling 
interval (VSI) and variable sample size (VSS) sampling designs for adaptive charts. Among all 
the advanced developments in the statistical process control (SPC) field, the CUSUM chart is 
an efficient alternative to the Shewhart chart when it comes to the detection of small to moder- 
ate shifts, see [1]. 

Both memory-type and memory-less control charts are used to monitor and control the 
spread or location parameters in any process. The sensitivity of memory-type charts have 
attracted the attention of many researchers and brought many more improvements in the SPC 
field and management of many organizations. Due to the efficiency of the CUSUM and 
EWMaA charts, many other control charts have been developed. Lucas and Crosier [2] incor- 
porated a fast initial response (FIR) feature to the classical CUSUM and compared its perfor- 
mance to that of the classical CUSUM chart. Zhao et al. [3] developed a dual CUSUM scheme 
which combines two CUSUM chart to detect a small shift in a process. Liand Wang [4] pro- 
posed an adaptive CUSUM Q chart also known as the self-starting approach. Shafqat et al. [5] 
introduced a linear prediction model based on the double EWMA control chart for non-nor- 
mal situations for monitoring the location parameter. Haq [6] and Над et al. [7, 8] proposed 
the use of EWMA control charts to monitor various processes, i.e. the mean and variance 
using ranked set sampling and the mean under the effect of measurement errors using ranked 
set sampling. Sanusi et al. [9] introduced a different approach for the EWMA-based chart 
when supplementary information about the main variable is not constant. Some different 
applications of attributes and variables control charts can be found in [10, 11]. 

The classical Shewhart control charts cannot be used efficiently when data are collected 
from uncertain and random situations, or when there are inconsistencies in the collected data. 
Khademi and Amirzadeh [12] reported that “fuzzy data exist ubiquitously in the modern 
manufacturing processes”; hence, in certain monitoring environments, ‘fuzzy model’ control 
charts are preferred when parameters or observations are uncertain, see for instance [13]. Sen- 
turk and Erginel [14] introduced the dispersion control chart by using the fuzzy logic, and 
Faraz et al. [15] proposed a control chart combining uncertainty and randomness assump- 
tions. Smarandache [16] provided an overview of the fuzzy logic and reported that the neutro- 
sophic logic deals efficiently with indeterminacy. In [17, 18], the neutrosophic logic is used to 
develop the idea of the neutrosophic statistics. Neutrosophic statistics is an extension of the 
classical statistics that is more efficient in analyzing data from ambiguous environment. For 
publications on neutrosophic statistics applied in various fields of knowledge, the readers are 
referred to [19-21]. For non-normal processes, control charts with a fuzzy approach are dis- 
cussed by [22]. Pereira et al. [23] designed a fuzzy chart for monitoring the blood components 
under the violation of the normality assumption. The most recent research works on control 
charts that use neutrosophic statistics concepts, can be found in [24-28]. 

At the best of the authors’ knowledge, there is no work on the CUSUM control chart under 
neutrosophic statistics. In this paper, the CUSUM X control chart using the neutrosophic sta- 
tistical method is proposed. Following the reports in [16, 17], it is presumed that the proposed 
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neutrosophic CUSUM (NCUSUM) X chart will provide very interesting properties (i.e. will be 
more efficient in detecting a small shift in the process and more informative, flexible, and 
effective in uncertain conditions). 

The rest of this paper is structured as follows: The second section gives a brief description 
of the classical CUSUM X chart and an introduction to the NCUSUM X control chart. In the 
third section, the performance of the proposed NCUSUM X chart is extensively investigated 
in terms of the neutrosophic run length properties. Moreover, the performance of the pro- 
posed chart is compared with the classical CUSUM X chart. Fourth and fifth sections provide 
illustrative applications of the new chart using simulated and real-life datasets, respectively. 
The conclusion and recommendations are given in the last section. 


The CUSUM control chart 


The classical CUSUM chart. The CUSUM X chart for monitoring small and moderate 
shifts in the mean was first introduced by [29]. The two-sided CUSUM chart is used to moni- 
tor upwards and downward shifts іп the process parameters using the charting statistics C} 
апа Сг, respectively, where i represents the sample number (or sampling time). Assume that X 
is a sequence of independent normally distributed observations given that the in-control mean 
Mo and variance т; thus, the charting statistics of the two-sided CUSUM X chart (i.e. С? and 
C;, with Су = Су = 0) are plotted versus a single control limit denoted by Н. These charting 
statistics are mathematically defined by 


ee (1) 





where X is the sample mean, ио is the target mean value of X and K is the reference value. The 
CUSUM chart depends on two singleton parameters, i.e. K and H, which need to be selected 
carefully because the sensitivity of the CUSUM chart depend on them, see pages 404-405 of 
Montgomery [1]. These parameters are chosen such that: 


К = ky/Var(X) and H = hy/Var(X), (2) 


where Var(X) = о2 /п, о, is the variance of X, k is selected to be half of the shift (i.e. К =°/,, 
with 6 is the shift in standard deviation units) and h is usually selected to be approximately five 
times the standard deviation units. Note that the process has moved upwards if C} > Н for 
any value of i; however, if C7 > Н, the process has moved downwards. 

In the next sub-section, the proposed NCUSUM X chart is introduced using the framework 
of the classical CUSUM X chart. 

The proposed neutrosophic CUSUM chart. In the uncertain environments having 
inconsistency data, the quality characteristic of interest (denoted by Ху) is not a deterministic 
value; hence, it falls within some interval, say [Хть Xvi], where the subscripts L and U refers to 
the lower and upper bound of those particular values; for instance, Ху; is a lower bound of Xy; 
and Хү; is the upper bound of Xy; (i.e. Xz;<Xn;<Xy;). The latter lower and upper notation 
holds for other expressions defined in this subsection. That is, the neutrosophic quality charac- 
teristic of interest Ху; Хм:Є[Хіь Хут], follows neutrosophic normal distribution with neutro- 
sophic mean fons иом [иг Hu] and neutrosophic variance 02; оўу Є [02,02], with и <иџ 
and a? < 02. Hence, the sample mean of the variable Xy; is given by Ху, Є [Ху Ху), 

Var(X pi) = бз / пу; ny Є [n ny] with np<ny, ом:Є[01 Oui] is the observed neutrosophic 
standard deviation оЁХ, with o,;<oy;. Thus, the plotting statistics of the proposed NCUSUM 
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X chart under uncertain environments is given by 


My: = max|0, Xy; — (Hoy + Ky) + Му) 


(3) 





М = тах[0, (Moy — Ky) — а Муу] 
where МУ € [M], М], Mx, € [М, М5], KwE [Ki Ku], НуЄ[Нү„ Ну] and pon is the target 
mean value. In addition, руу = HontO0on, Where Mon€ [Mor Hou] and ит [ит Hiu] represents 
the in-control and out-of-control neutrosophic mean, respectively; with 5 denoting the shift in 
the mean. Moreover, 


Ky = ky4/ Var(X y) and Hy = hy4/ Var(X y); ky € [ki kyl; hy € [hi hyl, (4) 


where ky and hy are the design parameters of the proposed NCUSUM X chart. The initial val- 
ues for the plotting statistics which are given in Eq (3) are equal to zero, i.e. М = My, = 0. 
The plotting statistics of the proposed chart are plotted versus the control limit Hy. The evalua- 
tion rule for the proposed chart is given as: for any value of i, if the value of Му, > Hy then the 
process of the mean is declared to have shifted upward; otherwise, if the value of My, > Hy 
then the process mean is declared to have moved downward. The values of ky and hy need to 
be selected with caution because the run length (RL) properties of the chart depend on them. 
Note that the RL is the number of charting statistics that are required before the first out-of- 
control signal is issued by a control chart. The manner in which the design parameters ky and 
hy are obtained is discussed in the next section. 


Run length performance 


Performance of the NCUSUM X chart. For some selected values of ky and ny, for exam- 
ple ky€{[0.20,0.25], [0.4,0.5]}, and ny€{[3,5],[10,20]}, the corresponding values of hy are 
found by running 10° simulations in R software for some given fixed nominal neutrosophic 
average RL (denoted Бу NARL») values of 300, 370, 400, 500. Tables 1-3 display the in-control 
and out-of-control neutrosophic average and standard deviation of the RL (denoted as NARL 
апа NSDRL) profiles of the proposed chart using different parameters, where NARLE[NARL,, 
NARLy] and NSDRLE[NSDRL,, NSDRLy], with NARL,<NARLy and NSDRL,<NSDRL y. In 
other words, Tables 1-3 give the ideal combinations of hy and ky when the nominal NARLo 
values are 300, 370, 400 and 500 for different values of ny. 

The following algorithm was applied to determine the values of ky and ny. 


Step-1: specify the values of NARLo and ny. 


Step-2: Determine the values of ky and ny where NARL>NARLp . Chose the values of ky and 
ny where NARL is very close or exact to NARL». 


Step-3: Determine the values of NARLs and NSDRLs for various values of 6. 


The performance of the proposed NCUSUM X control chart is evaluated in terms of differ- 
ent RL criteria, і.е. NARLs, NSDRL and the percentiles of the run length (PRL), і.е. the 25, 
50° and 75° percentiles which are denoted by P25 (і.е. P25€[P25;, P25y]), P50 (i.e. P50€ 
[Р50;, P50y]) and P75 (i.e. Р75Є[Р75;, P75y]). The following analysis would be noticed 
through Tables 1-6: 


1. Tables 1-5 show that there are inverse relationships between the values of out-of-control 
NARL (denoted as МАРКІ) and the values of ô. That is, the smallest value of NARL, is 
found at the largest values of 6, i.e. [1.00,1.00]. In other words, as 6 increase, the values of 
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Table 1. Тһе NARL and NSDRL values of the proposed NCUSUM X chart when ny€[3, 5], ky€[0.20, 0.25] for different nominal NARL. 


















































NARL 300 370 400 500 

hy [7.602,8.662] [8.099,9.299] [8.319,9.499] [9.039,10.299] 

ô NARLs NSDRLs NARLs NSDRLs NARLs NSDRLs NARLs NSDRLs 

0 [301.07,301.18] | [271.92,283.10] | [370.28,371.70] | [330.58,342.12] | [400.39,401.79] | [360.15,371.83] | [500.44501.54] | [441.99,447.12] 
0.25 [22.74,32.70] [12.50,19.00] [24.94,35.23] [13.38,20.73] [25.64,36.61] [13.68,20.68] [27.04,40.06] [14.01,22.71] 
0.5 [9.11,13.77] [3.38,5.35] [9.77,14.57] [3.32,5.38] [10.20,15.04] [3.49,5.47 [11.15,15.74] [3.61,5.58 
0.75 [5.00,8.59 [1.61,2.56] [6.36,9.33] [1.75,2.80] [6.54,9.52] [1.77,2.66 [6.96,10.14] [1.76,2.78 

1 [4.24,6.27 [1.06,1.60] [4.68,6.67] [1.09,1.69] [4.82,6.87] [1.06,1.71 5.20,7.36] [1.15,1.71 
1.5 [2.99,4.21 [0.58,0.87] [3.17,4.51] [0.57,0.87] [3.24,4.51] [0.61,0.89 3.44,4.86] [0.61,0.91 

2 [2.25,3.22 [0.44,0.57] [2.44,3.38] [0.50,0.60] [2.46,3.42] [0.51,0.59 2.67,3.70] [0.51,0.60 
2.5 [1.99,2.60 [0.33,0.51] [2.03,2.81] [0.29,0.56] [2.03,2.84] [0.20,0.56 2.14,3.03] [0.34,0.49 

3 [1.87,2.18 [0.15,0.38] [1.96,2.29] [0.19,0.45] [1.96,2.38] [0.17,0.46 1.99,2.60] [0.29,0.39 

4 [1.15,1.97 [0.16,0.35] [1.26,1.99] [0.08,0.43] [1.36,2.00] [0.09,0.44 1.62,2.01] [0.18,0.30 

5 [1.00,1.57 [0.00,0.29] [1.00,1.80] [0.03,0.39] [1.00,1.85] [0.05,0.40 1.02,1.96] [0.15,0.25 

6 [1.00,1.07 [0.00,0.16] [1.00,1.20] [0.00,0.20] 1.00,1.22] [0.01,0.22 1.00,1.55] [0.00,0.23 

7 [1.00,1.00 [0.00,0.03] [1.00,1.00] [0.00,0.07] [1.00,1.00] [0.00,0.09 1.00,1.03] [0.00,0.12 
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NARL, decrease until they approach a value of one, і.е. [1.00,1.00]. The latter implies that, 
for large values of 6, the control chart will give an out-of-control signal on the next sam- 
pling point. Similarly, the NSDRLs tend to zero for largest values of 6, і.е. [0.00,0.00], which 


means that the variability is reduced when 6 values are large. 


For a specified neutrosophic sample size, ny, there is an increasing trend in NARL, values 


as the neutrosophic value of ky increase at small shifts of б; however, there is an decreasing 
trend in NARL, values for large shifts of 5 (see Tables 1 and 2). For instance, when 
kn€[0.20,0.25], nnE[3,5] and 6 = 0.25 for a NARLo = 370, the NARL,€[24.94,35.23]; how- 
ever, when ky€[0.40,0.50], the NARL,€[28.53,41.29] using the same parameters. Note 
though, at ô = 3 for the same NARLy value, the NARL, is equal to [1.96,2.29] and [1.08,1.82] 


Table 2. Тһе NARL and NSDRL values of the proposed NCUSUM X chart when ny€[3, 5], ky€[0.4, 0.5] for different nominal NARL. 




















































































































NARI 300 370 400 500 
hy [4.585,5.421] [4.787,5.751] [4.887,5.894] [5.179,6.314] 
ô NARLs NSDRLs NARLs NSDRLs NARLs NSDRLs NARLs NSDRLs 
0 [301.20,301.83] | [292.12,299.37] | [370.88,371.32] | [354.69,355.21] | [400.82,401.85] | [370.48,374.30] | [501.32,501.62] | [450.81,452.90] 
0.25 [25.76,35.74] [20.24,27.67] [28.53,41.29] [22.35,32.16] [28.51,43.22] [21.76,33.36] [31.94,48.29] [25.25,37.55] 
0.5 [8.07,12.05] [4.07,6.65] [8.27,12.84] [4.07,6.84 [8.38,13.11] [2.98,6.43 [9.13,13.70] [4.53,6.49] 
0.75 4.61,6.67 [1.74,2.60] [4.70,7.09 [1.69,2.72 [4.87,7.20 [1.77,2.68 5.08,7.65] [1.78,2.91] 
1 3.27,4.79 [0.99,1.64] [3.35,5.06 [0.99,1.64 [3.41,5.15 [1.04,1.69 [3.55,5.46] [1.04,1.63] 
1.5 2.17,3.04 [0.47,0.81] [2.24,3.15 [0.52,0.81 [2.25,3.28 [0.52,0.79 [2.35,3.51] [0.55,0.86] 
2 1.74,2.29 [0.45,0.51] [1.81,2.43 [0.40,0.53 [1.85,2.45 [0.40,0.55 [1.91,2.61] [0.36,0.58] 
2.5 1.32,1.98 [0.32,0.46] [1.36,2.03 [0.32,0.48 [1.43,2.05 [0.30,0.32 [1.49,2.12] [0.30,0.37] 
3 1.04,1.74 [0.19,0.43] [1.08,1.82 [0.27,0.38 [1.08,1.86 [0.27,0.34 [1.15,1.93] [0.25,0.30] 
4 1.00,1.12 [0.00,0.33] [1.00,1.21 [0.00,0.20 [1.02,1.26 [0.00,0.30 [1.00,1.42] [0.09,0.19] 
5 1.00,1.00 [0.00,0.00] [1.00,1.00 [0.00,0.07 [1.00,1.01 [0.00,0.09 [1.00,1.02] [0.00,0.06] 
6 1.00,1.00 [0.00,0.00] [1.00,1.00 [0.00,0.00 [1.00,1.00 [0.00,0.00 [1.00,1.00] [0.00,0.00] 
7 1.00,1.00 [0.00,0.00] [1.00,1.00 [0.00,0.00 [1.00,1.00 [0.00,0.00 [1.00,1.00] [0.00,0.00] 
https://doi.org/10.1371/journal.pone.02461 85.1002 
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Table 3. The NARL and NSDRL values of the proposed NCUSUM X chart when ny€[10, 20], к< [0.4, 0.5] for different nominal NARL. 


















































NARLo 300 370 400 500 

hy [5.478,7.389] [4.804,5.730] [4.917,5.859] [5.224,6.172] 
ô NARLs NSDRLs NARLs NSDRLs NARLs NSDRLs NARLs NSDRLs 
0 [301.10,301.60] | [282.31,301.61] | [370.27,370.54] | [353.60,359.17] | [401.29,401.47] | [374.27,385.41] | [500.08,501.92] | [453.29,460.42] 

0.25 [12.32,13.55] [3.84,6.87] [8.43,14.49] [4.28,8.06 [8.50,14.96] [4.02,7.65 [9.09,15.13] [4.12,7.73] 

0.5 4.89,5.46 [0.97,1.80] [3.45,5.51 [1.06,1.94 [3.47,5.65 [1.01,1.88 [3.67,6.11] [1.08,2.07] 

0.75 3.16,3.29 [0.49,0.94] [2.23,3.50 [0.51,0.90] [2.27,3.62 [0.49,0.96 [2.39,3.74] [0.56,0.99] 
1 2.32,2.86 [0.45,0.61] [1.81,2.67 [0.41,0.61 [1.84,2.68 [0.40,0.64 [1.92,2.79] [0.33,0.65] 
15 1.63,1.87 [0.23,0.34] [1.08,1.95 [0.27,0.37 [1.09,1.94 [0.28,0.58 [1.17,2.00] [0.25,0.37] 
2 1.02,1.14 [0.00,0.27] [1.00,1.41 [0.03,0.25 [1.00,1.49 [0.00,0.50 [1.00,1.58] [0.00,0.49] 

2.5 1.00,1.00 [0.00,0.14] [1.00,1.06 [0.00,0.20 [1.00,1.04 [0.00,0.20 [1.00,1.09] [0.00,0.28] 
3 1.00,1.00 [0.00,0.00] [1.00,1.00 [0.00,0.03 [1.00,1.00 [0.00,0.04 [1.00,1.00] [0.00,0.03] 
4 1.00,1.00 [0.00,0.00] [1.00,1.00 [0.00,0.00 [1.00,1.00 [0.00,0.00 [1.00,1.00] [0.00,0.00] 
5 1.00,1.00 [0.00,0.00] [1.00,1.00 [0.00,0.00 [1.00,1.00 [0.00,0.00 [1.00,1.00] [0.00,0.00] 
6 1.00,1.00 [0.00,0.00] [1.00,1.00 [0.00,0.00 [1.00,1.00 [0.00,0.00 [1.00,1.00] [0.00,0.00] 
7 1.00,1.00 [0.00,0.00] [1.00,1.00 [0.00,0.00 [1.00,1.00 [0.00,0.00 [1.00,1.00] [0.00,0.00] 
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for ky equal to [0.20,0.25] and [0.40,0.50], respectively. Moreover, when design parameters 
are kept fixed, as ô increase, there is a decrease in the values of NARL}. 


From Tables 2 and 3, it can be noticed that the NARL, decreases rapidly when ny increase, 


with the restriction for the value of [k;, ky] being constant. For instance, when 6 = 0.25 and 
[Кт„ ky] = [0.40,0.50] for a nominal NARLy = 370, the NARL, is equal to [28.53,41.29] and 
[8.43,14.49] when ny is equal to [3,5] and [10,20], respectively. This indicates improving 
detection ability. 


The in-control neutrosophic RL distribution of the proposed chart is positively skewed as 


the in-control NARL are greater than the in-control P50. For instance, for the design 
parameters ky€[0.20,0.25] and ny€[3,5] corresponding to ô = 0, the in-control values of 


Table 4. The neutrosophic percentile points of the proposed NCUSUM X chart when ny€[3, 5], kv€[0.20, 0.25] for different nominal NARL. 












































NARLo 300 370 400 500 

ô P25 P50 P25 P75 P50 P75 P25 P50 

0 [98,104] [220,230] [111,123] [399,420] [273,277] [502,518] [119,132] [293,294] 
0.25 [14,19] [21,28] [15,21] [28,41] [22,30] [31,44] [16,22] [23,31] 
0.5 [7,10] [9,13] [7,11] [11,17] [9,14] [12,18] [8,11] [10,14] 
0.75 5,7 [6,8 [5,7] [7,10] [6,9] [7,11] [5,8] [6,9 

1 4,5 [4,6 [4,5] [5,7 [5,6] [5,8 [4,6] [5,7 
15 3,4 [3,4 [3,4] [3,5 [3,4] [3,5 [3,4] [3,4] 

2 2,3 [2,3 [2,3] [3,4 [2,3] [3,4 [2,3] [2,3 
2.5 2,2 [2,3 [2,3] [2,3 [2,3] [2,3 [2,3] [2,3 

3 2,2 [2,2 [2,2] [2,2 [2,2] [2,3 2,2] [2,2 

4 1,2 [1,2 [1,2] [1,2 [1,2] [2,2 [1,2] [1,2 

5 1,1 [1,2 [1,2] [12 [1,2] [1,2 [1,2] [1,2 

6 1,1 [1,1 [1,1] [1,1 [1,1] [1,1 [1,1] [1,1 

7 1,1 [1,1 [1,1] [1,1 [1,1] [1,1 [1,1] [1,1 





















































https://doi.org/10.1371/journal.pone.02461 85.1004 
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Table 5. The neutrosophic percentile points of the proposed NCUSUM X chart when ny€[3, 5], ку [0.4, 0.5] for different nominal NARL. 















































NARLo 300 370 400 500 
ô P25 P50 P25 P75 P50 P75 P25 P50 
0 [87,91] [201,202] [109,109] [405,428] [264,268] [504,528] [120,120] [285,293] 
0.25 [12,17] [19,27] [12,18] [34,47] [22,33] [37,54] [13,19] [22,33] 
0.5 [5,7 [7,10] [5,8] [10,15] [7,11] [10,16] [6,8] [7,12] 
0.75 [3,5 4,6 [3,5] [6,8] [4,7] [6,8] [4,5] 5,7 
1 [3,4 3,5 [3,4] [4,6 [3,5 [4,6] [3,4] 3,5 
1.5 [2,3 2,3 [2,3] [2,3 [2,3 [3,4] [2,3] 2,3 
2 [1,2 2,2 [2,2] [2,3 [2,2] [2,3] [2,2] 2,2 
2.5 [1,2 1,2 [1,2] [2,2 [1,2] [2,2] [1,2] 12 
3 [1,1 12 [1,2] [1,2 [1,2] [1,2] [1,2] 12 
4 [1,1 1,1 [1,1] [1,1 [1,1 [1,1] [1,1] 1,1 
5 [1,1 1,1 [1,1] [1,1 [1,1] [1,1] [1,1] 1,1 
6 [1,1 1,1 [1,1] [1,1 [1,1] [1,1] [1,1] 1,1 
7 [1,1 1,1 [1,1] [1,1 [1,1] [1,1] [1,1] 11 
























































https://doi.org/10.1371/journal.pone.0246185.t005 


Table 6. The neutrosophic percentile points of the proposed NCUSUM X chart when ny€[10, 20], к< [0.4, 0.5] for different nominal NARL. 


the P50€[273,277] are smaller than NARL, = 370 (see for instance Tables 1-6). The latter 


implies that the distribution of the NCUSUM X chart is positively skewed. 


5. Regardless of the NARLo value, the sensitivity of the NCUSUM X chart increase rapidly as 6 
increases in terms of the P25, P50 and P75 values, see Tables 4-6. 


Comparison of the proposed NCUSUM X and classical CUSUM X charts. The persis- 
tent need for more awareness of the proposed NCUSUM X chart, comprehensive comparison 
has been made with classical CUSUM X chart. To get valid deductions, the in-control ARLs of 
both charts are set at nominal levels of 370 and 500, and the performance is compared in terms 
of the out-of-control ARLs and SDRLs. The strategy of working in this section is parallel to 
[27], that is, a control chart under neutrosophic design is said to be capable if it has smaller 
NARL values compared to its counterpart. Thus, the efficiency of the proposed chart in terms 








































































































МАКЕ 300 370 400 500 
ô P25 P50 P25 P75 P50 P75 P25 P50 
0.00 [89,90] [201,205] [111,114] [408,408] [250,251] [511,547] [120,121] [275,276] 
0.25 [5,8 [7,12] [5,9] [10,17] [8,12] [10,18] [6,9] [8,13] 
0.5 [3,4 3,5 [3,4] [4,6] [3,5 [4,6] [3,4] 3,5 
0.75 [2,3 2,3 [2,3] [2,4 [2,3 [3,4] [2,3] 2,3 
1 [1,2 2,2 [2,2] [2,3 [2,3 [2,3] [2,2] 2,3 
1.5 [1,2 1,2 [1,2] [1,2] [1,2] [1,2] [1,2] 1,2 
2 [1,1 1,1 [1,1] [1,2] [1,1] [1,2] [1,1] 11 
2.5 {1,1 1,1 [1,1] [1,1] [1,1] [1,1] [1,1] 1;1 
3 {1,1 1,1 [1,1] [1,1 [1,1] [1,1] [1,1] 11 
4 {1,1 1,1 [1,1] [1,1] [1,1] [1,1] [1,1] 1,1 
5 [1,1 1,1 [1,1] [1,1 [1,1] [1,1] [1,1] 1,1 
6 [1,1 1,1 [1,1] [1,1] [1,1] [1,1] [1,1] 1,1 
7 [1,1 1,1 [1,1] [1,1 [1,1] [1,1] [1,1] 1,1 
https://doi.org/10.1371/journal.pone.0246185.t006 
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Table 7. Comparison of the NCUSUM and the classical CUSUM (CCUSUM) X charts when К = 0.5, kye[0.4, 0.5]. 















































CCUSUM chart: ARL = NCUSUM chart: NARLo = 370 CCUSUM chart: ARLo = NCUSUM chart: NARLy = 500 
370 500 
һ = 4.822 ny€[3, 5], һуЄ[4.887, 5.894] h=5.296 ny€[3, 5], hye [5.224, 6.172] 
ô ARLs SDRLs NARLs NSDRLs ARLs SDRLs NARLs NSDRLs 
0 370.11 354.62 [370.88,371.32] [354.69,355.21] 500.86 458.86 [501.32,501.62] [501.32,501.62] 
0.25 125.27 122.52 [28.53,41.29] [22.35,32.16] 162.99 150.86 [31.94,48.29] [48.29,31.94] 
0.5 35.37 28.00 [8.27,12.84] [4.07,6.84] 43.01 33.77 [9.13,13.70] [13.70,9.13] 
0.75 16.39 10.48 4.70,7.09 [1.69,2.72] 18.15 11.89 [5.08,7.65] 7.65,5.08 
1 9.98 5.25 3.35,5.06 [0.99,1.64] 11.46 6.37 [3.55,5.46] 5.46,3.55 
1.5 5.60 2.18 2.24,3.15 [0.52,0.81] 5.96 2.19 [2.35,3.51] 3.51,2.35 
2 3.85 1.23 1.81,2.43 [0.40,0.53] 419 1.27 [1.91,2.61] 2.61,1.91 
2.5 3.00 0.85 1.36,2.03 [0.32,0.48] 3.26 0.89 [1.49,2.12] 2.12,1.49 
3 2.48 0.63 1.08,1.82 [0.27,0.38] 2.72 0.71 [1.15,1.93] 1.93,1.15 
4 1.95 0.37 1.00,1.21 [0.00,0.20] 2.07 0.38 [1.00,1.42] 1.42,1.00 
5 1.62 0.48 1.00,1.00 [0.00,0.07] 1.80 0.40 [1.00,1.02] 1.02,1.00 
6 1.24 0.43 1.00,1.00 [0.00,0.00] 1.40 0.49 [1.00,1.00] 1.00,1.00 
7 1.04 0.19 1.00,1.00 [0.00,0.00] 111 0.31 [1.00,1.00] 1.00,1.00 












































https://doi.org/10.1371/journal.pone.02461 85.1007 


of the NARL values is compared to the ARL values of the classical CUSUM X chart. For 


instance, the classical CUSUM X chart with parameters k = 0.5 and h = 4.822 identifies a shift 
of size equal to 0.25 after 125.27 charting statistics for an АКТ» of 370; however, the proposed 
NCUSUM X chart with parameters ky€[0.4,0.5] and hy€[4.887,5.894] identifies a shift of size 
equal to 0.25 between the interval [28.53,41.29] samples for an МАКЕТ» of 370. In general, the 
comparison in Table 7 indicates that for small, moderate and large shifts in the process mean, 
the NARL and NSDRL are smaller in magnitude as compared to the ARL and SDRL of the clas- 
sical chart, at each corresponding shift value; and thus, the proposed NCUSUM X chart per- 
forms better than the classical CUSUM X chart. In summary, the neutrosophic and classical 
ARLs and SDRLs of the NCUSUM and CUSUM X charts indicate that the NCUSUM X chart 
has the higher shift detection ability as compared to the classical CUSUM X chart. 


Application of the NCUSUM X chart using simulated data 


In this section, the performance of the NCUSUM X chart is compared to one of the classical 
CUSUM X chart using simulated data. Fifty observations are generated using the neutrosophic 
normal distribution. The first twenty observations are generated by assuming that the process 
is statistically in-control and the next 30 observations are generated by assuming that the pro- 
cess has shifted with ô = 0.25. The simulated data along with ny€[3,5], к< [0.4,0.5], 
hy€[4.887,5.894], donE[1.00,1.00], and NARL,€[28.53,41.29], so the shift should be detected 
between the 28" and 41“ samples. The results from the simulated data are displayed in Figs 1 
and 2 for the proposed NCUSUM chart and the existing classical CUSUM X chart, 
respectively. 

From Figs | and 2, it can be seen that the proposed NCUSUM X chart detects the first shift 
at around the 30" sample in between the 28" and 41“ samples; however, the classical CUSUM 
X chart indicates that all the charting statistics are in-control. Hence, it can be deduced that, 
the proposed NCUSUM X chart has the ability to detect the shifts in the process mean earlier 
as compared to its classical counterpart. Therefore, the NCUSUM X chart is more efficient 
and performs better than the classical CUSUM X chart. 
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0 10 20 30 40 50 
Fig 1. The proposed NCUSUM X chart for simulated dataset. 
https://doi.org/10.1371/journal.pone.0246185.g001 





Real-data application of the proposed CUSUM X chart 


This section presents case-studies for the NCUSUM X chart that were applied to two types of 
data; the first example was dedicated to the case-study on the Pakistan state oil (PSO), while 
the second example is the case-study on the weather forecasting. The data were taken from the 
Stock Exchange and Meteorological Department of Pakistan. However, as mentioned in [28], 
all the observations and measurements of the current variables are not crisp numbers; they are 


CUSUM 
0 





О 10 20 зо 40 50 
Fig 2. The classical CUSUM X chart for simulated dataset. 





https://doi.org/10.1371/journal.pone.0246185.g002 





PLOS ONE | https://doi.org/10.1371/journal.pone.0246185 February 4, 2021 9/16 


PLOS ONE 


CUSUM control chart under uncertainty 





neutrosophically varied with some indeterminacy in representing intervals instead of singleton 
numbers. 

Petroleum industry case-study (PSO). The oil prices are strongly related with the 
stock prices, which is dependent on the principle of demand and there are some indeterminacy 
in its measurements. The dataset in Table 8 denotes the stock price in 8 different days of 
August of the last 10 years for the PSO. For this data, let ny = 8 and МАКТ = 370. The refer- 
ence value Ky = 0.50 y/./My (i.e. К = 0.5), and decision interval Hy = hoyy/./ny, where 
Oon€ [1.9221,2.239], ny = 8 and h = 4.804 are determined by using the procedure used to derive 
the values in Tables 1-3. Thus, the calculated decision interval and reference value given in Eq 
(4) is Hy€[3.265,3.803] and Ky€[0.3397.0.3958]. The calculated charting statistics of the 
NCUSUM chart (і.е. МУ and Мт) according to Eq (3) are provided in the last two columns of 
Table 8 and also plotted in Fig 3; however, those of the existing classical CUSUM X chart are 
plotted in Fig 4. 

From Fig 3, it can be noted that the PSO petroleum stock price are out-of-control in 2013 
to 2015 (і.е. higher prices); while, in 2011 and 2012 the petroleum oil prices are in-control 
(within reasonable range). However, after 2015, the prices are decreasing and are in-control. 
From Fig 4, it can be seen that in all the monitored years, the PSO petroleum prices are in-con- 
trol. From the comparison of both the charts, it can be seen that the proposed neutrosophic 
chart performs better and detects shifts earlier as compared to the classical counterpart. 

Pakistan meteorological case-study. For the Meteorological Department dataset of Paki- 
stan, the application of the proposed NCUSUM X chart is discussed. Since weather warning is 
a vital forecast because it is used for protecting people’s property and well-being as well as 
some socio-economic benefits. Thus, numerical weather experts use advanced mathematical 
formulas to predict the weather. Table 9 shows the weather forecast data information from 9 
different cities in Pakistan for 12 months. The abbreviations symbols for the cities name are: 
Lahore = LHR, Karachi = KHI, Multan = MUX, Faisalabad = FSD, Islamabad = ISB, 

Peshawar = PEW, Hyderabad = HDD, Skardu = KDU and Gilgit = GIL. 

Let the sample size and the in-control ARL values be ny = 12 and МАКІ, = 370, while the 
reference and decision interval values are Hy€[1.4109,2.6919] апа Ky€[0.1472,0.2809] when 
h = 4.791, k = 0.5, and оомуЄ[1.0197,1.9464]. So, the calculated NCUSUM charting statistics 
(ie. Mx, and M,,) are given in the last two columns of Table 9. The charting statistics of the 
proposed NCUSUM X chart are shown in Fig 5 and the ones of the classical CUSUM X chart 
are shown in Fig 6. 

From Fig 5, the proposed NCUSUM X chart is found to be more efficient and gives better 
results compared to its classical counterpart in Fig 6; as it indicates downwards shift in the 
weather forecast and shows that the weather is very cold in Skardu and Gilgit cities. Note 
though, the classical CUSUM X chart reveals that all cities’ weather are within normal range 
and does not detect the shift at any point in the data. Thus, in comparison of both the control 
charts for monitoring weather forecasting, we can say that the NCUSUM X chart provides bet- 
ter prediction especially when meteorologists are not sure about the temperature values to be 
used for weather prediction. Therefore, the neutrosophic CUSUM X chart is recommended 
for weather prediction under the uncertainty environments. 


Conclusion and future recommendations 


In this paper, a new NCUSUM X control chart is proposed. The neutrosophic run length per- 
formance of the proposed chart is computed using Monte Carlo simulations. The application 
of the proposed chart is illustrated using simulated data and two real-life examples using the 
PSO petroleum industry dataset and the weather forecasting dataset from nine cities in 
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Fig 3. Real-life data from PSO petroleum oil company’s stock price for the proposed NCUSUM X chart with 
Hy€[3.265, 3.803] and Ky€[0.3397, 0.3958]. 


https://doi.org/10.1371/journal.pone.0246185.g003 
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Fig 4. Real-life data from PSO petroleum company’s stock price for the classical CUSUM X chart with h = 4.804 and К = 0.5. 
https://doi.org/10.1371/journal.pone.0246185.g004 
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Fig 5. The proposed NCUSUM X chart for the weather condition in Pakistan with Hy€[1.4109, 2.6919] and Ky€[0.1472, 
0.2809]. 


https://doi.org/10.1371/journal.pone.0246185.g005 
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Fig 6. The classical CUSUM X chart for the weather data with h = 4.791 and К = 0.5. 
https://doi.org/10.1371/journal.pone.0246185.g006 
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Pakistan. The results show that the NCUSUM X chart is better suited for use than the classical 
CUSUM X chart either when the data are inaccurate or when the operators do not have a com- 
plete dataset or for uncertain environments situations. The comparison study shows that the 
proposed neutrosophic chart is efficient in detecting small to large shifts in the process param- 
eter than the classical chart. Note that the NCUSUM X chart is functional when the variable of 
interest adopts the neutrosophic normal distribution. The proposed neutrosophic chart is also 
recommended for monitoring the process mean shifts for the medical instruments, automo- 
biles, aerospace and drinking water industries. 

Further future works can be dedicated to the NCUSUM X chart that can be designed for 
non-normal distributions, as well as for the multiple dependent state sampling and repetitive 
sampling approaches. 


Acknowledgments 


The authors are deeply thankful to the editor and reviewers for their valuable suggestions to 
improve the quality of the original manuscript. 


Author Contributions 

Data curation: Ambreen Shafqat. 

Funding acquisition: Mohammed Albassam. 
Software: Ambreen Shafqat. 

Supervision: Muhammad Aslam. 
Visualization: Mohammed Albassam. 
Writing – original draft: Muhammad Aslam. 


Writing – review & editing: Muhammad Aslam, Jean-Claude Malela-Majika, Sandile С. 
Shongwe. 


References 


1. Montgomery DC, Introduction to Statistical Quality Control, 6th Edition, John Wiley & Sons, Hoboken, 
New Jersey; 2009. 


2. Lucas JM, Saccucci MS. Exponentially weighted moving average control schemes: Properties and 
enhancements. Technometrics. 1990; 32(1): 1—12. https://doi.org/10.1080/00401706.1990.10484583 


3. Zhao YI, Tsung Е, Wang F. Dual CUSUM control schemes for detecting a range of mean shifts. ПЕ 
Transactions. 2005; 37(11): 1047—1057. https://doi.org/10.1080/074081 70500232321. 


4. 12, Wang Z. Adaptive CUSUM of the О chart. International Journal of Production Research. 2010; 48 
(5): 1287—1301. https://doi.org/10.1080/00207540802484937. 


5. Shafqat A, Huang Z, Aslam М, and Nawaz MS. A Nonparametric repetitive sampling DEWMA control 
chart based on linear prediction. IEEE Access; 2020; 8: 74977-74990. https://doi.org/10.1109/ 
ACCESS.2020.2989132. 


6. Haq A. An improved mean deviation exponentially weighted moving average control chart to monitor 
process dispersion under ranked set sampling. Journal of Statistical Computation and Simulation. 
2014; 84(9), 2011—2024 https://doi.org/10.1080/00949655.2013.780059 


7. Haq A, Brown J, and Moltchanova Е. New exponentially weighted moving average control charts for 
monitoring process mean and process dispersion. Quality and Reliability Engineering International. 
2015; 31(5): 877-901. https://doi.org/10.1002/qre.1646. 


8. Haq A, Brown J, Moltchanova Е, Al-Omari Al. Effect of measurement error on exponentially weighted 
moving average control charts under ranked set sampling schemes. Journal of Statistical Computation 
and Simulation. 2015; 85(6): 1224—1246. https://doi.org/10.1080/00949655.2013.873040. 





PLOS ONE | https://doi.org/10.1371/journal.pone.0246185 February 4, 2021 15/16 


PLOS ONE 


CUSUM control chart under uncertainty 





10. 


11. 


12. 


13. 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


21. 


22. 


23. 


24. 


25. 


26. 


27. 


28. 


29. 


Sanusi RA, Riaz М, Adegoke NA, Xie М. Ап ЕММА monitoring scheme with a single auxiliary variable 
for industrial processes. Computers & Industrial Engineering. 2017; 114: 1-10. https://doi.org/10.1016/ 
j.cie.2017.10.001. 


Shafqat A, Huang Z, Aslam М. Design of X-bar control chart based on Inverse Rayleigh Distribution 
under repetitive group sampling. Ain Shams Engineering Journal. 2020; https://doi.org/10.1016/j.asej. 
2020.06.001. 


Shafqat A., Hussain J., Al-Nasser A. D., and Aslam M. Attribute control chart for some popular distribu- 
tions. Communications in Statistics—Theory and Methods. 2018; 47(8):1978—1988. https://doi.org/10. 
1080/03610926.2017.1335414. 


Khademi М, Amirzadeh V. Fuzzy rules for fuzzy x and r control charts. Iranian Journal of Fuzzy Sys- 
tems., 2014; 11(5): 55—66. https://doi.org/10.22111/IJFS.2014.1722. 


Sabegh MHZ, Mirzazadeh A, Salehian S, Weber GW. A literature review on the fuzzy control chart: 
classifications & analysis. International Journal of Supply and Operations Management. 2014; 1(2): 
167—189. https://doi.org/10.22034/2014.2.03. 


Senturk S, Erginel N. Development of fuzzy Х—Ё and X —$ control charts using a-cuts. Information Sci- 
ences. 2009; 179(10):1542—1551. https://doi.org/10.1016/j.ins.2008.09.022. 


Faraz A, Kazemzadeh RB, Moghadam MB, Bazdar A. Constructing a fuzzy Shewhart control chart for 
variables when uncertainty and randomness are combined. Quality & Quantity, 2010; 44: 905-914. 
https://doi.org/10.1007/s11135-009-9244-9. 


Smarandache F. Neutrosophic Logic—A generalization of the intuitionistic fuzzy logic. SSRN Electronic 
Journal, 2016. https://dx.doi.org/10.2139/ssrn.2721587. 


Smarandache F. Introduction to Neutrosophic Statistics, Sitech & Education Publishing, Craiova, 
Romania & Ohio, USA; 2014. 


Khalid HE, Essa AK. Introduction to Neutrosophic Statistics, the translated Arabic version, Pons Pub- 
lishing House, Brussels, Belgium, 2020. 


Moran JCS, Chuga JFE, Collaguazo WMA, Gudiño CWM. Neutrosophic statistics applied to the analy- 
sis of socially responsible participation in the community. Neutrosophic Sets and Systems, 2019; 
26:19-28. https://doi.org/10.5281/zenodo.3244232. 


Maldonado PAC, Martinez YP, Valverde GSE, Erazo JDI. Neutrosophic statistics methods applied to 
demonstrate the extra-contractual liability of the state from the Administrative Organic Code. Neutro- 
sophic Sets and Systems, 2019; 26:29-34. https://doi.org/10.5281/zenodo.3244262. 


Carballido RM, Paronyan H, Matos MA, Molina ALS. Neutrosofic statistics applied to demonstrate the 
importance of humanistic and higher education components in students of legal careers. Neutrosophic 
Sets and Systems, 2019; 26:174—180. https://doi.org/10.5281/zenodo.3244848 


Panthong C, Pongpullponsak A. Non-normality and the fuzzy theory for variable parameters X control 
charts. Thai Journal of Mathematics, 2016; 14(1):203—213. 


Pereira P, Seghatchian J, Caldeira B, Xavier S, de Sousa G. Statistical methods to the control of the 
production of blood components: Principles and control charts for variables. Transfusion and Apheresis 
Science, 2018; 57(1):132—142. https://doi.org/10.1016/j.transci.2018.02.022 PMID: 29526479 


Aslam M, Bantan RAR, Khan N. Design of a new attribute control chart under neutrosophic statistics. 
International Journal of Fuzzy Systems, 2019; 21:433-440. https://doi.org/10.1007/s40815-018-0577- 
1. 


Aslam M, Rao GS, Shafqat A, Ahmad L, Sherwani RAK. Monitoring circuit boards products in the pres- 
ence of indeterminacy. Measurement, 2021; 168:108404. https://doi.org/10.1016/j.measurement. 
2020.108404. 


Khan Z, Gulistan M, Hashim R, Yaqoob N, Chammam W. Design of S-control chart for neutrosophic 
data: An application to manufacturing industry. Journal of Intelligent & Fuzzy Systems. 2020; 38(4): 
4743—4751. https://doi.org/10.3233/JIFS-191439. 


Aslam M. Attribute Control Chart Using the Repetitive Sampling under Neutrosophic System. IEEE 
Access, 2019; 7:15367—15374. https://doi.org/10.1109/ACCESS.2019.2895162. 


Aslam M. Control chart for variance using repetitive sampling under neutrosophic statistical interval sys- 
tem. IEEE Access, 2019; 7: 25253-25262. https://doi.org/10.1109/ACCESS.2019.2899020. 


Page ES. Continuous inspection schemes. Biometrika, 1954; 41(1—2): 100—115. https://doi.org/10. 
1093/biomet/41.1-2.100. 





PLOS ONE | https://doi.org/10.1371/journal.pone.0246185 February 4, 2021 16/16 


